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Abstract 

A Hamiltonian approach to the solution of the Vlasov-Poisson equations has 
been developed. Based on a nonlinear canonical transformation, the rapidly os- 
cillating terms in the original Hamiltonian are transformed away, yielding a new 
Hamiltonian that contains slowly varying terms only. The formalism has been 
applied to the dynamics of an intense beam propagating through a periodic fo- 
cusing lattice, and to the coherent beam-beam interaction. A stationary solution 
to the transformed Vlasov equation has been obtained. 
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1. INTRODUCTION 



The evolution of charged particle beams in accelerators and storage rings can often be 
described by the Vlasov-Maxwell equations^'^. At high energies the discrete-particle collision 
term^ comprises a small correction to the dynamics and can be neglected. Radiation effects 
at sufficiently high energies for leptons can be a significant feature of the dynamics, and 
should normally be included in the model under consideration. 

The Vlasov-Maxwell equations constitute a considerable simplification in the description 
of charged particle beam propagation. Nonetheless there are only a few cases that are 
tractable analytically^'^. Therefore, it is of the utmost importance to develop a systematic 
perturbation approach, able to provide satisfactory results in a wide variety of cases of 
physical interest. 

Particle beams are subject to external forces that are often rapidly oscillating, such as 
quadrupolar focusing forces, RF fields, etc. In addition, the collective self-field excitations 
can be rapidly oscillating as well. A typical example is a coUiding-beam storage ring device, 
where the evolution of each beam is strongly affected by the electromagnetic force produced 
by the counter- propagating beam"^^^. The beam-beam kick each beam experiences is local- 
ized only in a small region around the interaction point, and is periodic with a period of one 
turn. 

In this and other important applications, one is primarily interested in the long-time 
behavior of the beam, thus discarding the fast processes on time scales of order the period 
of the rapid oscillations. To extract the relevant information, an efficient method of aver- 
aging is developed in Sec. 2. Unlike the standard canonical perturbation technique^'^, the 
approach used here is carried out in a "mixed" phase space (old coordinates and new canon- 
ical momenta)^, which is simpler and more efficient in a computational sense. It should 
be pointed out that the formalism developed here is strictly speaking non-canonical and in 
general does not provide complete elimination of fast oscillating terms in the transformed 
Vlasov equation in the mixed variables. Nevertheless, such an elimination can be performed 
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in the new Hamiltonian in the mixed variables. Moreover, if the distribution function is 
assumed to be an arbitrary function of the new time-independent Hamiltonian, it is a sta- 
tionary solution of the nonlinear Vlasov equation in the mixed variables. The canonical 
perturbation method developed in this paper is further applied to intense beam propagation 
in a periodic focusing structure (Sec. 3) and to the coherent beam-beam interaction (Sees. 
4 and 5). A coupled set of nonlinear integral equations for the equilibrium beam densities 
is derived. 

To summarize, the effectiveness of the Hamiltonian formalism developed in the present 
paper is demonstrated in two particular examples. In the first example discussed in Sec. 
3, the short-scale dynamics is contained in the external focusing force acting on the beam, 
while an essential feature of the coherent beam-beam interaction treated in Sees. 4 and 5 is 
the relatively fast variation of the collective interaction between the colliding beams. The 
simplicity in applying the Hamiltonian averaging technique is embedded in the use of mixed 
canonical variables. Instead of expanding the generating function and the new Hamiltonian 
in terms of the new canonical coordinates and momenta*^'^ one has to simply solve the 
Hamilton-Jacobi equations order by order. It should be emphasized that the mixed variable 
Hamiltonian formalism can be used to derive amplitude equations, describing processes of 
formation of patterns and coherent structures in a number of plasma and beam systems in 
which collective processes are important. 



2. THE HAMILTONIAN FORMALISM 

We consider a A^-dimensional dynamical system, described by the canonical conjugate pair 
of vector variables (q, p) with components 

q= (?i,?2,---,gjv), 

P= (Pi,P2,---,Piv)- (1) 
The Vlasov equation for the distribution function /(q, p; t) can be expressed as 
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f + [/,^kp = 0, (2) 
where 

^ aF^G _ dFdG 

^ ' J^'P dqi dpi dpi dqi ^ ' 

is the Poisson bracket, i?(q, p; t) is the Hamiltonian of the system, and summation over 
repeated indices is imphed. Next we define a canonical transformation via the generating 
function of the second type according to 

5 = 5(q,P;i), (4) 
and assume that the Hessian matrix 

of the generating function 5'(q, P; t) is non-degenerate, i.e., 

det {Uij) ^ 0. (6) 

This imphes that the inverse matrix T-Lj^ exists. The new canonical variables (Q, P) are 
defined by the canonical transformation as 

We also introduce the distribution function defined in terms of the new canonical coor- 
dinates (Q, P) and the mixed pair of canonical variables (q, P) according to 

/o(Q, P; t) = /(q(Q, P; t), p(Q, P; t); t), (8) 



Fo(q,P;0 = /(q,p(q,P;i);i). (9) 

In particular, in Eq. (8) the old canonical variables are expressed in terms of the new 
ones, which is ensured by the implicit function theorem, provided the relation (6) holds. As 



far as the function Fo(q, P;i) is concerned, we simply replace the old momentum p by its 
counterpart taken from the first of Eqs. (7). Because 

we can express the Poisson bracket in terms of the mixed variables in the form 

(dF dG dF dG\ 
^^^^ 

Differentiation of Eq. (7) with respect to time t, keeping the old variables (q, p) fixed, yields 



dqidt dqidPj \ dt 



qp 



'dQi\ d^S d^S (dP, 



+ ^^^ ' (13) 



\dt dPidt dPidP. \ dt , 

\ ' qp J \ ' qp 

or 

Our goal is to express the Vlasov equation (2) in terms of the mixed variables (q, P). 
Taking into account the identities 

dqj ~ dqjdPi ~ ^ dQj ~ ' ^ ^ 

dfo dFo dfo d^S 



dPi dPi dQj dPidPj ' 

we obtain 

(OA ^dfo df^fdQA dfofm] 
[dt)^^ dt dQ^dt)^^ dpAdt)^^ 

dt ^'[dqidtdPj dPjdtdqi) 



(17) 
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dt 



+ 



dS_' 
'dt 



(18) 



q,P 



Furthermore, using the relation 
where 

n{q,P;t)^H{ci,W,S;t), 

we express the Vlasov equation in terms of the mixed variables according to 

dt 
where 



(19) 



(20) 



(21) 
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/C(q, P;t)^— + Hid, V,S; t) 



(22) 



is the new Hamiltonian. 

For the distribution function /o(Q,P;i), depending on the new canonical variables, we 
clearly obtain 

dfo 



dt 



+ [/o,^]q,p = 0. 



(23) 



where the new Hamiltonian /C is a function of the new canonical pair (Q, P), such that 
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)C{VpS,P;t) = — + H{q,V,S;t), 



(24) 



and the Poisson bracket entering Eq. (23) has the same form as Eq. (3), expressed in the 
new canonical variables. 



3. PROPAGATION OF AN INTENSE BEAM THROUGH A PERIODIC FOCUS- 
ING LATTICE 



As a first application of the Hamiltonian formalism, we consider the propagation of a con- 
tinuous beam through a periodic focusing lattice in a circular ring with radius R. Particle 
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motion is accomplished in two degrees of freedom in a plane transverse to the design orbit. 
The model equations consist of the nonhnear Vlasov-Poisson equations^ 

i + [/'^U = 0> (25) 

Vl7j;^-ing^-injd'pf{q,p;9), (26) 
where 

^(q, p; ^) = f {pI +pi) + jR (^-^' + ^^^') + ^^(q; ^) (27) 

is the normalized Hamiltonian, and q = {x,z). The transverse canonical momenta p = 
{Px-iPz) entering the Hamiltonian (27) are dimensionless variables which represent the actual 
transverse momenta of the moving particle scaled by the longitudinal momentum of the 
synchronous particle^'^. 

In addition, R is the mean radius of the accelerator and ^0 is a normalized potential 
related to the actual electric potential tp according to 

where N is the total number of particles in the beam, is the particle charge, and Eq is the 
electric susceptibility of vacuum. Moreover, the parameter A is defined by 

A - ^, 

where fig — Vg/c is the relative velocity of the synchronous particle, 7^ = (1 — Pi) is the 
Lorentz factor, and 

is the classical radius of a beam particle with charge and rest mass m^. The coefficients 
Gx,z{d) determining the focusing strength in both transverse directions are periodic functions 
of ^ 



G,,,{9 + e)^G,,,{9), (31) 
with period ©. 

Following the procedure outhned in the preceding section we transform Eqs. (25) - (27) 
according to 

(32) 

— +ei/(q,Vq5;e) = /C(q,P), (33) 

Vl^iP = -47r J d^PFoiq, P) det (VqVp5), (34) 

where e is formally a small parameter proportional to the focusing field strength, which will 
be set equal to unity at the end of the calculation. The next step is to expand the quantities 
S, K, and in a power series in e according to 

So + €Si + e^S2 + e^-Sg + . . . , (35) 
/C = /Co + e/Ci + e^/Ca + e^/Cg + . . . , (36) 

V' = ^0 + e^i + £^^-2 + £^3 + ■ ■ ■ ■ (37) 

We now substitute the expansions (35) - (37) into Eqs. (33) and (34) and obtain perturbation 
equations that can be solved order by order. 

The lowest order solution is evident and has the form 

5-0 = q • P, /Co = 0, (38) 

V^V^o = -47Tjd^PFo{q,P). (39) 

First order 0(e): Taking into account the already obtained lowest order solutions (38) 
and (39), the Hamilton- Jacobi equation (33) to first order in e can be expressed as 



Imposing the condition that the first order Hamiltonian fCi be equal to 



we obtain immediately 



(40) 



(41) 



(42) 



V'l = 0. 

Here we have introduced the notation 

'o+e 

e 



Oo+e 



Gx,M- j dT[GxM-G: 



(43) 



(44) 



Note that since the focusing coefficients are periodic functions of 9 they can be expanded 
in a Fourier series 



GxM = E Ggexp(ml]^), 

n=— oo 

where 



(45) 



1 



e 



G^::l = ^ j d9GxMexp{-tnn9), 





(46) 



and Q = 2n/Q. Therefore for the quantities Gx,z and Gx,zi9) expressed in terms of the 
Fourier amplitudes, we obtain 



<^.,.W = -^E^exp(znO^). 



Gx,. = Gg), GxAO) = - o E -r^ exp {inn9). (47) 

Second order 0{e^): To this order, the Hamilton- Jacobi equation (33) takes the form 



dS2 
d9 



- (xPxGx + zP.G,) =/C2(q,P). 



(48) 



It is straightforward to solve Eq. (48), yielding the obvious result 



/C2(q,P) = 0. 



(49) 



For the second order potential ■02 we obtain the equation 
= -47r(§, + I d^PFo(q, P), 



(50) 



or, making use of (39), 



G,{9) + G,{9) 



V'o(q)- 



(51) 



In Eqs. (49) - (51), Gx,z{9) denotes apphcation of the integral operation in Eq. (44) to 
Gx,z{^): i.e., 



G,M = [ dTG,,z{T), 



(52) 



because G^^z — 0. 

Third order 0{e^): To third order in e, the Hamilton- Jacobi equation (33) can be written 



as 



^ + Rl^p'^d^ + p^dz) + ^{cy + cy) + a(4 + dz^A = iCsiq, P). 

The third-order Hamiltonian IC^ is given by the expression 



(53) 



(54) 



Equation (53) can be easily solved for the third-order generating function S3. The result is 

1 



-R[ P^Gx + PzGz 



— i^Glx' + Glz']-X[Gx + Gz]^|^o. 



(55) 



For the third-order electric potential ip^ we obtain simply 



^3 = 0. 



(56) 



Fourth order 0(e'^): To the fourth order in the expansion parameter e the Hamilton- 
Jacobi equation (33) can be expressed as 
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^-xP. {gI + G^G,) -zP, [gI + G.G^i -Xr(g, 



dx 



dz 



/C4(q,P). 
(57) 



The obvious condition to impose is that the fourth-order Hamiltonian /C4 be equal to 



^4(q, P) — —xPxGxGx — zPzGzGz- 



(58) 



With Eq. (58) in hand, it is straightforward to solve the fourth-order Hamilton- Jacobi 
equation (57) for S^. We obtain 

-54 = xPx {&x + Gxb)j + zPz {dl + Gzb)j + \R {dx + [P.^ + Pz^^ . (59) 
For the fourth-order electric potential ■04, we obtain the Poisson equation 



GxGz + G2 + + GxGx + GzGz + \R \Gx + J ^0 



(60) 



Fi/t/i order O(e^): In fifth order, we are interested in the Hamiltonian /C5. Omitting 
algebraic details we find 

/C5(q, P) = f + gA + 1 (g.S^x^ + Gz^lz' 



+ A 



9a; 



9^ 



(61) 



In concluding this section, we collect terms up to fifth order in e in the new Hamiltonian 
K, — 1Cq + e/Ci -|- e^lC2 + ■ ■ ■ and set e = 1. This gives 



^(q,P)= E {^Pl + BuuPu + ^u']+\U^) 



u={x,z) 



c 



+ A 



dx 



dz 



where the coefficients Au, Bu and C„ are defined by the expressions 



(62) 



(63) 
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(64) 



(65) 



The Hamiltonian (62), neglecting the contribution from the self-field V'O) describes the un- 
perturbed betatron oscillations in both horizontal and vertical directions. 

It is useful to compute the unperturbed betatron tunes v^^z in terms of averages over 
the focusing field-strengths. For a Hamiltonian system governed by a quadratic form in 
the canonical variables of the type in Eq. (62), it is well-known that the characteristic 
frequencies yx,z can be expressed as 



y"^ - AC - 



{u — x,z). 



(66) 



Keeping terms up to sixth order in the perturbation parameter e, we obtain 



(67) 



In terms of Fourier amphtudes of the focusing coefficients, the Eq. (67) can be expressed as 



^2 ^ ^(0) ^j^y\ - 



2^ 

fi4 



m,n=l 



Re(Gi™)*GWGi™-")) 
mn{m — n)^ 



-2 E 

m:ri=l 



Re(G(,™)GWG'(;"+")*) 
mn{m -\- nf 



+ 



4e6 - \G^r^ Gi") 



E 



m,n=l 



(68) 



For purposes of illustration, we consider a simple FODO lattice with equal focusing and 
defocusing strengths -\-G and — G, and period ©. We also assume that the longitudinal 
dimensions ^/ of the focusing and defocusing lenses are equal; the longitudinal dimensions 
Od of the corresponding drift spaces are assumed to be equal as well. Moreover, 



2{ei + Od) = e. 



(69) 



For simplicity we consider the horizontal degree of freedom only (the vertical one can be 
treated in analogous manner). The Fourier amplitudes of the focusing coefficients are 
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G(2n+i) = ^^^^^{cxp [-.(2n + l)nef] - 1}, Gi^n) ^ 0^ (70) 

where n = 0, 1, 2, . . .. To second order in e, we obtain for the horizontal betatron tune 



.2 



m= 



TT^ ^1 (2m - 1)4 e 
In the hmit of infinitely thin lenses, Of — > 0, Eq. (71) reduces to the well-known expression 

= (72) 
where use of the identity 

°° 1 TT^ 

has been made. 

It is evident from Eqs. (68) and (71), that the Hamiltonian averaging technique developed 
here represents a powerful formalism for evaluating the betatron tunes in terms of averages 
over the focusing field strength. 

4. COHERENT BEAM-BEAM INTERACTION 

As a second application of the Hamiltonian formalism developed in Sec. 2, we study here 
the evolution of two counter-propagating beams, nonlinearly coupled by the electromagnetic 
interaction between the beams at collision. For simplicity, we consider one-dimensional 
motion in the vertical (g) direction, described by the nonlinear Vlasov-Poisson equations 

de 



+ [/,,i7fe]=0, (74) 



where 



' 47rldpfs.k{q,p;e), (75) 



Hk = -f{p' + q')+hSMyk{q;0) (76) 
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is the Hamiltonian. Here Xk is the beam-beam couphng parameter, defined according to^^ 

_ RreN3_kP*kg 1 + /5fcO/5(3-fc)0 2RreN3-kP*k, 
"^k — p ^2 ~ F • V ' ' J 

7fc0-t>(3-fc)a; Pfco lk0'l^{3-k)x 

Moreover, {k — 1,2) labels the beam, fk{q,P',d) is the distribution function, 9 is the az- 
imuthal angle, and i/k is the betatron frequency in vertical direction. In addition, R is the 
mean machine radius, r^. is the classical electron radius, Ni^2 is the total number of par- 
ticles in either beam, Vk{q;9) is the normalized beam-beam potential, is the vertical 
beta-function at the interaction point, and Lj^j. is the horizontal dimension of the beam 
ribbon^^. 

Our goal is to determine a canonical transformation such that the new Hamiltonian is 
time-independent. As a consequence, the stationary solution of the Vlasov equation (21) is 
expressed as a function of the new Hamiltonian. Following the procedure outlined in the 
Sec. 2 we transform Eqs. (74) - (76) according to 

F^'\)Ck\=0, (78) 
^ + eH,(^q,^;e^^}Ck(q,n (79) 

where e is again a formal small parameter, which will be set equal to unity at the end of the 
calculation. 

The next step is to expand the quantities Sk, IC^ and in a power series in e, analogous 
to Eqs. (35) - (37), according to 

S,^qP + 6G« + e^Gi^) + e^Gf + . . . , (81) 

/C, = 6/C« + + + . . . , (82) 

Vk^V, + eVt'^ + 6V,(^) + e'vt'^ + . . . , (83) 
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where 

^^AnldPFt'\<l,n (84) 

Substitution of the above expansions (81) - (83) into Eqs. (79) and (80) yields perturbation 
equations that can be solved successively order by order. The results to third order in e are 
briefly summarized below. 
First Order. 0(e) 

}Ci\^,P)-'j(P' + <l') + ^m, (85) 

Gi\q, P; 9) = '-^Mq) E - exp {tnO), (86) 

y«(g;^) = 0. (87) 
Second Order. O(e^) 

/Cf(?,P) = 0, (88) 

(g, P; 9) = -^PV;'(g) E A ^xp (m^), (89) 

(?; ^) = -^V^f (?) E exp (zn^), (90) 



where 

Third Order. 0{e^) In third order we are interested in the new Hamiltonian, which is of 
the form 



4nVl'{q)JdPFt'\q,P). (91) 



^'i\Q,P) = '-^m[V^{q) - 2V^''{q)\, (92) 
where (^{z) is Riemann's zeta-function^^ 
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OO -1 

n=l "' 



(93) 



The effectiveness of the Hamiltonian formahsm developed in the present paper has been 
demonstrated in two particular examples. In the first example discussed in the previous 
section, the short-scale dynamics is contained in the external focusing force acting on the 
beam, while an essential feature of the coherent beam-beam interaction treated above is the 
relatively fast variation of the collective interaction between the two colliding beams. The 
simplicity in applying the Hamiltonian averaging technique is embedded in the use of mixed 
canonical variables. Instead of expanding the generating function and the new Hamiltonian 
in terms of the new canonical coordinates and momenta^'® one has to simply solve the 
Hamilton- Jacobi equations order by order. It should be pointed out that the mixed variable 
Hamiltonian formalism can be used to derive amplitude equations, describing processes of 
formation of patterns and coherent structures in a number of plasma and beam systems in 
which collective processes are important. 

5. THE EQUILIBRIUM DISTRIBUTION FUNCTION 

Since the new Hamiltonian ICk is time-independent (by construction), the equilibrium dis- 
tribution function Fq'^^ [see Eq. (78)] is a function of the new Hamiltonian 



Integrating Eq. (94) over P we obtain a nonlinear integral equation of the Haissinski type 
for the equilibrium beam density profile Qq'^ 



(94) 



where 




(95) 




(96) 



where 
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Z{u,v) = sgn(M)sgn(w) — 2|m|(5(w). 



(97) 
(98) 
(99) 



Here sgn(2;) is the well-known sign-function. 

Let us further specify the function Qk{]^k) and assume that it is given by the thermal 
equilibrium distribution^'^'^'^^ 



(100) 



where Hk is a normalization constant, defined according to 



exp 



(101) 



and £fe is the unnormalized beam emittance. The second term in the Hamiltonian (97) can 
be transformed according to 



/| (3-fc) 

q-q\Qo 



{q') = g - (qs-k) + 2 / d?'(?' - q)Qt''\q') 



= q - (gs-fc 
where 



) +2 JdqiqiQ^o~^\qi + q), 



(102) 



CO 

(g.")= jdqqyt\q). 



(103) 



Expanding the beam density *^''(q'i + q) occurring in the integral in Eq. (102) in a Taylor 
series and integrating by parts, we obtain 

oo Q(n) 
n=2 

where 



^ oo Mn) 

/ dq'\q - q'\Qt'\q') = {qt\) ' + (1 - '^M-k)q + 2 E 

t/ Q 77/. 



(104) 
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in 



dqQQo\q), 



(105) 



in) _ d-'gi'\q) 



q=0 



Qqn-2 

u 

Substituting (100) and (104) into Eq. (96) we obtain 



(106) 



V 



(107) 



Taking into account that 



,(n) 



/27r£:A 



X. ->{ Vk \("-2)/2 



where Hn{z) is the Hermite polynomiaP^ of order n, we obtain 



i7„_2(0) + O(Afe), (108) 



Qi\q) = AfkJ^^ exp [/ife(g)], 



(109) 



where 



and 



Vkq^ XkBz-kq 27rAfce3_fcA/'3-fe 



2£i 



27r£feA4 



q^{a3-kq) + 



(110) 



(111) 



Here, $(^) is the error function^^. 

In order to determine the normahzation constant(s), j\4, we utihze the method of Laplace 
to take the integral of the beam density Qq^^ (g) over q. The first step consists in finding the 
extremum value(s) q^^^ of the function(s) hk{q)- They satisfy the (two) equation(s) 



, XkBs-k , 27rAjk£3_fej\/'3_A,.^^^ (e) 



+ 



+ 



£fe ^k £k^3-k 

These are evidently maxima, since 



aa-kqk 



0. 



;ii2) 



^ ^ Sk Skiys-k 



(113) 
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Integrating the beam density (109) over g, we obtain^^ 



27rM; 



exp 



(114) 



Equation (114) represents two transcendental equations for determining the normahzation 
constants A4- For the beam centroid and the beam size, i.e., the first and the second 
moments of the beam density (109), we obtain 

2A4 



exp 



hk{qt') 



(115) 



(Ik + ■ ■■/ (.\\ \/ — exp 



3/2 



exp 



hk{q^:^) 



(116) 



In order to proceed further, we assume that the beam-beam couphng parameter 
is small, and expand the equilibrium beam density go'\q) in a perturbation series in 
according to 



Qo \q) = ^00 (?) + hgoiil) + 



(117) 



where 

^00 (?) 
and 



k\/T^ 



exp (^-a\(f) , 



(118) 



^01 (?) = -— "( Os-fcg + — 



g$(a3_fcg) + 



^o?(?)- 



(119) 



The main goal in what follows is to determine the normalization constant (s) A4o- To do so 
we integrate Eq. (117) over q. As a result of simple algebraic manipulations, we obtain 

7rA4 TTi/TrAfc 



2 3 3 ^al + ai_feA4AA3_fe = 1. 

"ifc ^kU^k^^S-k 



Introducing the new unknowns 



(120) 
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= ^, (121) 
we can write the two equations for determining A4i,2 as 
1 = Ml - biMiM2, 

1 = M2 - b2MiM2, 
where 



(122) 



^ ^ ^ (123) 

£iV7r 0102 £2V^ 01^2 



Prom Eq. (122), as a result of simple algebraic manipulations we obtain the quadratic 
equation 

b2Ml - (62 -bi + l)Mi + 1 = (124) 
for J\4i, and the equation 

biM2 = b2Mi + 61-62 (125) 

for determining AI2 once Aii is known. Equation (124) has one real double root if and only 
if the discriminant 

P = (62 - 61 + if - 462 (126) 
is equal to zero. This gives 

62 = (y^^^il)'. (127) 
Since the scaled normalization constants M.1^2 should be positive we choose 

^61-1)'. (128) 
Thus we obtain 
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Ml = ^ = , M2 = (129) 



|v^-i|' Vh' 

To conclude this section we note that in the case of X> 7^ we have two solutions for 
either M., i.e., 

262 26i 
Note also that the discriminant "D is invariant (does not change) under permutation of bi 
and 62- In other words, four different physically realizable situations are possible for a wide 
range of parameters 

< 62 < 1 + (131) 

The inequality in Eq. (131) has been obtained under the condition that both solutions in 
Eq. (130) are positive. This case corresponds to the so-called "flip-flop" state^^ of the two 
colliding beams, which is a bifurcated state better to be avoided. 



6. CONCLUSIONS 

We have developed a systematic canonical perturbation approach that removes rapidly os- 
cillating terms in Hamiltonians of quite general form. The essential feature of this approach 
is the use of mixed canonical variables. For this purpose the Vlasov-Poisson equations are 
transformed to mixed canonical variables, and an appropriate perturbation scheme is chosen 
to obtain the equilibrium phase space density. It is worthwhile to note that the perturbation 
expansion outlined in the preceding section can be carried out to arbitrary order, although 
higher-order calculations become very tedious. 

In conclusion, it is evident from the present analysis that the Hamiltonian averaging 
technique developed here represents a powerful formalism with applications ranging from 
beam propagation through a periodic focusing lattice (Sec. 3) to coherent beam-beam 
interaction (Sees. 4 and 5). For example, in the application to the coherent beam-beam 
interaction, the rapidly oscillating terms due to the periodic beam-beam kicks have been 
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averaged away, leading to a new time-independent Hamiltonian (Sec. 4). Furthermore, 
the equihbrium distribution functions have been obtained as a general function of the new 
Hamiltonian, and coupled set of integral equations for the beam densities has been obtained 
(Sec. 5). An intriguing feature of the analysis in Sec. 5 is the derivation of a condition for 
existence of the so-called "flip-flop" state^^ of the two coUiding beams, which is a bifurcated 
state better to be avoided in experimental applications. 

We reiterate that the formalism developed here is strictly speaking non-canonical and in 
general does not provide complete elimination of fast oscillating terms in the transformed 
Vlasov equation in the mixed variables. Nevertheless, such an elimination can be performed 
in the new Hamiltonian in the mixed variables. Moreover, if the distribution function is 
assumed to be an arbitrary function of the new time-independent Hamiltonian, it is a sta- 
tionary solution of the nonlinear Vlasov equation in the mixed variables. 

Finally, we reiterate that the mixed variable Hamiltonian formalism developed in the 
present analysis can be used to derive amplitude equations, describing processes of formation 
of patterns and coherent structures in a number of plasma and beam systems in which 
collective processes are important. 
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